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Abstract. - The Kolmogorov stochasticity parameter is shown to act as a tool to detect point 
sources in the cosmic microwave background (CMB) radiation temperature maps. Kolmogorov 
CMB map constructed for the WMAP's 7-year datasets reveals tiny structures which in part coin- 
cide with point radio and Fermi/LAT gamma-ray sources. In the first application of this method, 
we identified several sources not present in the then available OFGL Fermi catalog. Subsequently 
they were confirmed in the more recent and more complete 1FGL catalog, thus strengthening the 
evidence for the power of this methodology. 
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Introduction. — A number of point sources 
are identified in the temperature maps of cos- 
mic microwave background (CMB), including 
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ciency of the search of point sources is particularly impor- 
tant for obtaining a pure cosmological signal of high pre- 
cision at large multipoles, for the cross correlation studies 
with galaxy surveys, on baryonic oscillations, e.g. [2] [3]. 

We now involve for this aim a new descriptor, the Kol- 
mogorov stochasticity parameter which describes the de- 
gree of randomness of a given sequence 0] [5] . It has been 
already been applied to the CMB temperature datasets 
[B]. The resulting Kolmogorov CMB maps showed the 
potential of this approach in the separation of the Galac- 
tic disk from the cosmological signal, in the identification 
of anomalies such as the non-Guassian Cold Spot [7] [5]. 
Certain spots and regions noticed in the K-map can be 
compared with the studies with other descriptors [9]. 

Here we constructed and analyzed Kolmogorov's map 
for WMAP's 7-year temperature datgQ and obtained a 
list of 12 regions of a degree scale with anomalously high 
value of Kolmogorov's function (see below). Among those 



lr rhis study initially was based on the WMAP's 5-year data and 
7-year data were involved later upon their release: the results are 
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regions, 6 coincide with sources given in the catalog of the 
WMAP's point sources, where they are identified as ra- 
dio galaxies, while all 12 regions coincide with gamma-ray 
sources discovered by Fermi satellite's Large Area Tele- 
scope (LAT) [TU]; 4 among these regions in both catalogs 
are overlapping mutually. 

In the first version of this study (arXiv:1002.2155v), 
when only OFGL catalog of Fermi/LAT was available 
|llj . there were 3 regions which have practically identi- 
cal values both of the Kolmogorov's function and of the 
CMB temperature distribution but were absent in either 
of mentioned catalogs, i.e. of the radio and Fermi gamma 
sources. One might thus expect that they also are still 
unidentified sample of point sources. Indeed, the release 
of Fermi/LAT 1FGL catalog confirmed this prediction, i.e. 
all 12 CMB structures are identified in that catalog, thus 
showing the predictive power of the approach. 

The detected sources include blazars/flat spectrum ra- 
dio quasars (FSRQs) and active galaxy [TU], which are dis- 
tinguished by their multi-wavelength emission, from radio 
to gamma rays, and due to their superluminous nature are 
efficient in probing the distant Universe. 

The Kolmogorov analysis of the CMB map. — 

The Kolmogorov statistic and the stochasticity parameter 
[5 [Sj is defined for a sequence {-X^i, X2-> * • * 3 3£ n J- repre- 
sented in increasing order, so that the distribution func- 
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tion is F(x) = P{X < x} . Together with the empirical 
distribution function 

r o, x <x i; 

F n (x) = l k/n, X k < x < X k+1 ,k = 1,2,..., n- 1; 
[ 1, X n <x, 

(1) 

the stochasticity parameter is defined as A„ = 
y/nsup x \F n (x) — F (x)\. For any continuous function F 
lim„^ oc P{A„ < A} = $(A) , where 

$(A) = ^ (- x ) fe e_2fc2A2 » $ (°) = 0, A > , (2) 

fe=— oo 

the convergence is uniform and $ is independent onF ||]. 

The parameter A n and the Kolmogorov's statistic de- 
fine the degree of the randomness of the sequence, which 
can be tested via numerical experiments [12| . Namely, the 
stochasticity parameter gives a quantitative measure of 
the presence of correlations in a sequence of numbers inde- 
pendent of the underlying probability distribution within 
the interval of the probable values of A = 0.4-1.8. This cri- 
terion was applied to sequences following from the number 
theory [5] [T3], and, as mentioned, to a physical problem, 
i.e. the CMB maps representing distribution of the tem- 
perature values over the set of pixels. The latter resulted 
in the distribution of Kolmogorov's function $ over the sky 
which can reflect the properties of the large-scale matter 
inhomogeneities in the Universe i.e. on the distribution of 
the voids |BJ. 

The Kolmogorov's CMB maps were obtained for the 
Gaussian function for F and for WMAP several w-band 
datasets in [T3] [T5] , where the procedure of averaging over 
the pixel sets in Healpix [TB| representation is described, 
zy-band has the best resolution of all bands and is the least 
contaminated by the Galactic synchrotron radiation which 
is particularly important for us, as some of the studied re- 
gions are of not high Galactic altitude. In the distribution 
of $, it appears that the regions with its values higher 
than 0.999 form a sample of 12 regions. Those 12 regions 
were found via two procedures: the 3> sky was scanned 
with a step 0.5° with averaging of $ within the regions of 
1.5° radii, and with a step 0.2° and averaging within 1.0°. 
The radius of the circles are determined by the available 
angular resolution of WMAP maps, i.e. at smaller regions 
the number of pixels becomes too small, while for larger 
regions the effect of the point source non-Gaussianity is 
smeared. The study of simulated maps in P3] [H3 had 
revealed differences in the Kolmogorov's parameter in the 
WMAP's and generated pure Gaussian maps, as the role 
of the WMAP's noise in the behavior and value of that 
parameter has been studied. The numerical simulations 
were performed also via generation of a given signal or 
several signals with the subsequent comparison of their $- 
maps with the observed ones, in order to reveal the con- 
tribution of a given component. These studies had shown 



that Kolmogorov's parameter does reflect certain statisti- 
cal properties of the CMB maps which enable the study of 
non-Gaussianities. Direct quantitative comparison of the 
source detecting power of this approach and of other meth- 
ods, however, is not much informative, since Kolmogorov's 
parameter is reflecting different type of correlations. 

The results of both scans revealed the same 12 regions 
at |6| > 10° as shown in Fig. 1. Fig. 2 contains an 
example for the temperature and $ behavior in one of 
those regions, No. 04, of radius 5°. In Fig. 3 we represent 
the CMB temperature behavior in the 12 regions denoted 
when radio (R), Fermi/LAT (G) point source counterparts 
are available. All $ sources possess (5-function type peak 
of temperature. Table 1 contains the coordinates of the 
revealed regions, the coordinates of the identified radio 
sources according to the WMAP's point source catalog, 
and Fermi/LAT gamma sources, along with the accuracy 
of identifications in arc minutes. 

Conclusion. — We used the Kolmogorov's parame- 
ter to extract point sources in the CMB maps. The idea 
is that the radiation of different origin can have different 
statistic and this can be revealed by the Kolmogorov's dis- 
tribution. Numerical experiments for systems with built 
a priori different statistic showed that K-parametcr is in- 
deed sensitive to the randomness properties of systems 
|12j . Namely, the regions which might be of lower signif- 
icance by conventional methods and hence are absent in 
WMAP's catalog, can be easily outlined by the present 
approach. On the other hand, the increase in the angu- 
lar resolution of the CMB temperature maps will enable 
to apply this method for smaller scales, i.e. to obtain the 
averaged Kolmogorov function for large enough number of 
pixels, and hence to increase its source detecting power. 

The analysis performed for WMAP's 7-year w-band 
maps aimed the extraction of around 1 degree scale regions 
(namely, when averaged within 1.0° and 1.5° circles at 0.2° 
and 0.5° step scans, respectively), which were outlined in 
the K-map due to their high value of the Kolmogorov's 
function When compared with the catalog of the point 
sources of WMAP's corresponding map, 6 among those 
^-regions were identified with radio galaxies, while all 12 
regions were identified with Fermi/LAT gamma-ray 1FGL 
sources [TU|. The 3 detected anomal regions were absent 
in the previous 0FGL Fermi/LAT catalog [TT|. However 
as concluded then, since the CMB regions have practi- 
cally identical distributions of the function $ as well as 
of the CMB temperature, the latter regions can also be 
point sources of close nature. Present study confirmed 
that prediction and the interest for dedicated studies of 
CMB map's those regions in radio (see e.g. [17]) and other 
bands. 

At the available angular resolution of WMAP the 
method enables to identify at high confidence level only 
a sample of 12 regions, which is certainly smaller than the 
number of sources identified by WMAP's mask. Again, 
the reason, as described above, is in the available reso- 
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Fig. 1: The location of the 12 high <E> regions in the sky, i.e. those outside the Galactic disk with |6| > 10°. 
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Fig. 2: The CMB temperature map region of 5° radius centered on No. 04 (a), and its $ distribution (b). 



lution of the WMAP's maps. However, even at present 
resolution the method identifies sources not noticed by 
WMAP's procedure, which appear to be Fermi/LAT 
gamma sources [10]. Most of those sources appear to be 
blazar/quasars and active galaxy, i.e. multi-wavelength 
and supcrluminous objects, and we showed another possi- 
bility of their identification based on totally different idea, 
the statistic of the signal. The absence of some of the 
detected sources in the initial OFGL catalog and their ap- 
pearance in the next, 1FGL catalog, shows the predictive 
power of the approach. 

The Kolmogorov parameter thus enables to detect point 
sources in the CMB maps not noticed by other criteria. 
The applicability of the method will increase at higher 
resolution maps expected from the PLANCK mission. 

REFERENCES 

[1] Jarosik N., Bennett C.L., et al, arXiv:1001.4744, 
(2010) 



[2] Sawangwit U. et al, Mon.Not.Roy.Astr.Soc, 402 (2010) 
2228 

[3] Sawangwit U. et al, arXiv:0912.0511, (2009) 

[4] Kolmogorov A.N., G.ht.Ital.Attuan,, 4 (1933) 83 

[5] Arnold V., Nonlinearity, 21 (2008) T109 

[6] Gurzadyan V.G. & Kocharyan A. A, A&A,, 492 (2008) 
L33; ,40,4, 493 (2009) L61 

[7] Cruz M., Martinez-Gonzalez E. and Vielva P., to be 
published in the "Highlights of Spanish Astrophysics V", 
Proceedings of the VIII Scientific Meeting of the Spanish 
Astronomical Society (SEA), Santander, July 7-11, 2008; 
\arXiv: 090 1.7986 1 (2009) 

[8] Gurzadyan V.G., Starobinsky A. A., et al., A&A, 490 
(2008) 929 

[9] Rossmanith G., Raeth C, et al., MNRAS, 399 (2009) 
1921 

[10] Fermi-LAT Collaboration, arXiv.1002.2280, (2010) 
[11] Abdo A. A., Ackermann A., et al, ApJS, 188 (2010) 
405 

[12] Ghahramanyan T., et al, Mod. Phys. Lett. A, 24 (2009) 
1187 

[13] Arnold V.I., ICTP/2008/001, Trieste, ((2008)) 



p-3 



V.G.Gurzadyan et al. 




61 62 63 64 65 66 67 68 -23 -22 -21 -20 -19 -18 -17 -16 




53 54 55 56 57 58 59 60 61 -16 -15 -14 -13 -12 -11 -10 -9 




13 14 15 16 17 18 19 20 -23 -22 -21 -20 -19 -13 -17 -16 




Decl. (degree) Decl. (degree) 



Fig. 3: The distribution of the CMB temperature in the 12 regions (with smoothed, blue, line) ; (R) and (G) denote those with 
radio and Fermi/LAT gamma-source identifications, respectively. 
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Table 1. Coordinates of the 12 regions in WMAP's CMB map outlined by Kolmogorov's function "1>, along with the coordinates 
of the identified radio sources, with the accuracy of identification in arcminutes dA and source designations, and the same for 
Fermi/LAT gamma-ray sources. 
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